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ABSTRACT 

We present an empirical assessment of the use of broadband optical colours as age indicators 
for unresolved extragalactic clusters and investigate stochastic sampling effects on integrated 
colours. We use the integrated properties of Galactic open clusters as models for unres olved 
extragalactic clusters. The population synthesis code Starburst99 dLeitherer et al.l[T999 ) and 
four optical colours were used to estimate how well we can recover the ages of 62 well-studied 
Galactic open clusters with published ages. We provide a method for estimating the ages of 
unresolved clusters and for reliably determining the uncertainties in the age estimates. Our 
results support earlier conclusions based on comparisons to synthetic clusters, namely the 
(U— B) colour is critical to the estimation of the ages of star forming regions. We compare 
the observed optical colours with those obtained from Starburst99 using the published ages 
and get good agreement. The scatter in the (B— V)obseri;eci— (B— V),„ode/ is larger for lower 
luminosity clusters, perhaps due to stochastic effects. 

Key words: open clusters and associations: general — galaxies: star clusters — galaxies: 
stellar content — methods: data analysis 



1 INTRODUCTION 

Thousands of luminous young star clusters have been dis- 
covered in external galaxie s by the Hubble Spa ce Telescope 
(HST) (e.g., iHoltzman et al. 1992. 1996; Whitmore et al. 1993'; 
Whitmore & Schvyeized 1 1995. ; .Meurer et al., .1995.; Johnson et al.i 



! 19991 ; [Benedict et al.l l20oS~lKeel & Bornd I2OO3I; iHancock et all 
20031 ; IWhitmord |2003| ; IWeistrop et al] |2004 and references 



therein). Population synthesis analyses show that these clusters are 
sometimes very young, with ages of a few to a few tens of Myr. 
Such age information can be extremely valuable, because it pro- 
vides clues to cluster formation processes as well as cluster de- 
struction mechanisms. A radial gradient in the ages of clusters 
within a galaxy can tell us about gas inflow driven by a bar or 
an interaction. For interacting galaxies, comparison of cluster ages 
with dynamical models can provide information about star form- 
ing mechanisms. For example, gas compression along the Arp 107 
tidal arm may hav e caused a gradient of ages along the arm (e.g., 
ISmith et al.ir2005 j) , while in M51, a burst of cluster formation oc- 
curr ed at around the time of the latest passage of the compan- 
ion teastian et alj2005l ; lLee, Chandar. & Whitmorel2005l) . Cluster 
synthesis studies have led to the suggestion that there is a high 'in- 
fant mortality' in star clus t ers, with m any dissolving within 10 Myr 
(e.g., iBastian et alj|2005l ; [pall, Chandar, & Whitmoreii2005i) . Star 
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clusters are themselves sometimes clustered into complexes with 
characteristic radii of ~1 kpc (e.g., Zhang, Fall, & Whitmore 200l|; 
iLarse nl2004l) , and the more massive clusters tend to be located near 
the centre of these complexes, suggesting cluster merging. Accu- 
rate age dating of young clusters associated with 'ultraluminous 
X-ray sources' (ULXs) can help distinguish between stellar-mass 
and intermediate-mass (10 0—1000 M[-n) black holes for the origin 
of the X-ray emission (e.g., ISmithet"al]|2005bl) . 

Stellar population synthesis models have evolved dra- 
matically in the last several decades, since the pioneering 
works where data on globular clusters and individual gi- 
ant stars were used to fill gaps in t he evolu t ionary tracks 
to make early colour mode ls (e.g. iTinslevI 1 19681 1 19721 ; 
IStruck-Marcell & TinslevI Il978h . Recently, several grou ps have 
introd uced ev olutionary synthesis co des, e.g., JBruzua L & Chariot! 
( I99 3I) (B&C),lFritze-v . Alvensleben & Gerhard IT9941) (GALEV) 
Fioc & Rocca-Volmerange, (,1997) (PEGASE), and lLeitherer et al.l 
( 199^ (Starburst99). These new models vary in terms of their in- 
put physics, stellar spectral libraries and extinction laws. In gen- 
eral, however, there is good agreement among these models (e .g., 
ICharlot. Worthev. & Bressanlll996l ; IVazauez & Leitherej2005l) . 

Unfortunately, however, at the present time, it is not clear 
how accurate age estimates based on these models are for unre- 
solved extragalactic clusters. In most cases spectra of the targets 
are not available, so the amount of extinction and metal abun- 
dance are unknown. In these cases, one typically compares some 
combination of UV, optical, and/or near-IR broadband colours to 
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model cluster spectral energy distributions (SED) to simultane- 
ously estimate age and extinction assuming some metallici t y (e.g . 
Pasauali, de G rijs. & Gallaghe 3 I2OO3I : iHancock et alJl200l bool 
Smit h et al. 200^1 

There are several parameters that can affect the colour of a 
cluster other than age. Reddening plays a very important role, as 
does the chemical composition. However, in some age ranges and 
colours, reddening and chemical composition are degenerate with 
age. Furthermore, for low mass clusters (< 10''' M0), the observed 
integrated colours can be affected by stochas tic sampling of the 
initial mass fu nction (IMF) (see for exampl e, iCerviiio & Luridianal 
I2OO6 , 2004; Ce rviiio & Valls-Gabau J2OO3I and references therein). 
For example, the random addition of a small number of high mass 
stars will affect the integrated colour of a cluster. This suggests that 
the observed integrated colours of low mass clusters may not be 
good indicators of age when compared to the integrated colours of 
a stellar population model with a fully sampled IMF. 

It is not clear how well a particular set of colours can predict 
the age of a cluster. Are some colour sets better suited than others? 
What uncertainties can be expected because of the choice of colours 
in the comparison? What uncertainties can be expected because of 
the assumptions made in generating the model SEDs? 

Several authors have investigated the use of broadband colours 
as age indicators by co mparison of model colours to the colours 
of synthetic clusters (e.g. lOil de Paz & Madordl 20021 : 1 Anders et"ai] 
l2004l : lde Griis et aEll 20051) . We use an alternative method by com- 
paring model colours to the integrated colours of resolved Galactic 
open clusters (OCs). We compare the published integrated colours 
of well-studied OCs to a set of population synthesis models. These 
clusters have published ages previously determined by the turn-off 
of the zero-age main sequence on the H-R diagram. This work par- 
allels that o f Pessev et al, (2008) , who do a similar analysis of Mag- 
ellanic Cloud clusters using both optical and near-IR data. 

The present paper is organized as follows. In §2 we describe 
our sample of OCs. We describe the population synthesis models 
used in this study in §3 and the data analysis in §4. Predicting the 
ages of the OCs, predicting the amount of extinction and the effects 
of metallicity, stochastic sampling effects and cluster dissolution 
are discussed in §5. Finally, we summarize in §6. 



2 THE DATA SAMPLE 

We started with the set of Galactic open clusters from the 
WEBD^lWeb Base Donnees Amas) database operat ed at the In- 
stitute for Astronomy of the University of Vienna ( lMermillio3 
1 1995! ). The WEBDA database contains 379 OCs, with well- 
determined integrated (U-B), (B-V), (V-R), an d (V-I) colours 
and colour excesses, as found by several authors dLata et alj|2002l: 
Battinelli, Brandimarti, & Capuzzo-Dolcettal 1994 : Pandev et al 



19891 : ISpassova & Beavll 19851 : ISagar. Joshi. & Sinvhallll983l : lGrav 



19651) . We then culled the WEBDA sample t o include only the 



sample of well-studied OCs in lPaunzen & Net opil (2006), who es- 
tablished a list of 72 open clusters with the most accurate known 
parameters to serve as a standard table for testing isochrones and 
stellar models. The age uncertainties in lPaunzen & Netopill ^20061) 
were determined by measuring the standard deviation of all the 
published ages in the literature for each of the OCs. 

Not all of the 4 integrated colours were determined for each 
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of the OCs in the|p aunzen & Netopill Jiooel) standard set. To model 
observations with unknown dust extinction, we reversed the extinc- 
tion corrections to the published colours using the published values 
of E(B— V) and the conversions from E(B— V) to t he oth er colour 
excesses. We used the conversions in iLata et al.l ( |2002|) . namely 
E(U-B)=0.72 E(B-V)+0.05 E(B-V)^ E(V-R)=0.6 E(B-V), 
and E(V-I)=1.25 E(B-V). 

From the standard set we created sub-samples of OCs for each 
of the 4 optical colours and 10 different combinations of colours. 
When multiple integrated colour measurements were available, we 
adopted the most recently determined values. Our final data sets 
include OCs with mean age uncertainties of 19% and ages ranging 
from 8 Myr to 8.8 Gyr. The metallicities associated with this sample 
ranges from [Fe/H] ~l/7 solar to ^--^2.3 solar. 

Unfortunately, the WEBDA database does not give the un- 
certa inties on the total colour s for individual OCs. According to 
Saga r. Joshi. & Sinvhall ( Il983l) , the maximum uncertainty in the 
published integrated colours for the WEBDA data set is ±0.2 mag. 
One of the sources of uncertainty listed is the error in the redden- 
ing. Because we have used the published colour excesses to reverse 
the extinction corrections, we can neglect the extinction uncertainty 
in the colours. Removing this from the total uncertainty, assuming 
it was originally added in quadrature, the maximum uncertainty in 
colour is ~0.I4. We assume that all the measured colours have this 
maximum uncertainty. This assumption further allows us to make 
fair comparisons of both the accuracy and precision afforded by 
each colour in age estimation. 



3 THE MODEL CLUSTERS 

We used a set of ev olutionary synthesis models from the Star- 
burst99 (SB99) code dUeitherer et ai] |l999'). We used the new v5.1 
code, which in cludes the Padova asy mptotic giant branch (AGB) 
stellar models jVazguez" & LeithereJ l2005). The new version ac- 
counts for all stellar phases that contribute to the integrated light 
of a stellar population with arbitrary age from extreme UV to NIR. 
Strictly speaking, the Geneva tracks are more appropriate for mod- 
eling young clusters, less than 10 Myr, when O stars are present. 
Most of our sample OCs have ages greater than 100 Myr so for 
simplicity we only consider models with the Padova tracks. 

Our SB99 model spectral energy distributions (SEDs) were 
generated assuming a Kroupa initial mass function (IMF) (favors 
high mass stars) ( iKroupai .2002) with exponents of 1.3 and 2.3 
and mass ranges from 0.1 — 0.5 M© and 0.5 — 100 Mq re- 
spectively. We have also assumed instantaneous (single burst) star 
formation, and solar abundances. It has been demonstrated that 
adopting di fferent forms of the I MF has a minor impact on opti- 
cal colours jMacArthur et alj|2004) . so we have not explored var- 
ious IMFs. We also generated models with abundances less than 
(0.2 x) and greater than (2.5 x) solar The model SEDs were kA- 
dened from 0.0 m a g to 2.0 mag in 0.02 mag increments using 
the lCardelh. Clayton. & Mathisi ( 1 1 989) reddening law. Finally, the 
model SEDs were convolved with the Johnson and Rron-Cousins 
UBVRI filter bandpasses and the broadband optical colours were 
determined. The models span a range of ages from 1 Myr to 20 
Gyr. From 1 Myr to 1 Gyr we used a step size of 1 Myr; from 1.1 
to 20 Gyr, the step size was 100 Myr. Because the models will be 
compared to the integrated colours of resolved stellar populations, 
only the stellar contributions were included in the SEDs. Nebular 
emission lines have been shown to be important in the first 10^ yr 
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(U-B) 

Figure 1. (B— V) vs. (U— B) colour-colour diagram. The black curves are 
our SB99 models with solar abundances and varying E(B— V). From the 
bottom to the top, E(B— V)=0.0 (thick soUd curve), 0.5 (dashed curve), 1.0 
(dashed curve), and 1.5 (dotted curve) mags. The filled squares are the ob- 
served colours of the standard set of open clusters in the WEBDA sample. 
For each curve, age increases to the right from 1 Myr to 20 Gyr. Note that 
we have assumed the maximum uncertainty in colour for each OC. 

d Anders & Fritze-v. Alvenslebenll2003h and should be considered 
when studying unresolved stellar populations (see §5.3). 

Figure 1 plots the (B— V) vs. (U— B) colour-colour diagram. 
The black curves are our SB99 models with solar abundances and 
varying E(B-V). From the bottom to the top, E(B-V)=0.0, 0.5, 
1.0, and 1.5 mags. The filled squares are the observed colours of 
the standard set of open clusters in the WEBDA sample. For each 
curve, age increases to the right from 1 Myr to 20 Gyr. Note that 
we have assumed the maximum uncertainty in colour for each OC. 



4 DATA ANALYSIS 

This study is not intended as a test of the SB99 model's ability 
to reproduce the colours o f young clusters. This has a lready been 
demonstrated (see for e.g.. lVazauez & Leithere r 2005). We intend 
to test how well various colours and initial model assumptions can 
recover the ages of well studied Galactic open clusters. In addition 
to using the observed colours to estimate ages, we compared the 
published colours with the model colours calculated assuming the 
published ages to look for evidence of stochastic sampling effects. 

To determine the model predicted ages of the OCs we 
compared the observed integrated colours to each of the red- 
dened (and unreddened) SB99 model colours, for a single as- 
sumed metallicity. We used a mi r umization calculatio n (e.g., 
IPasauah. de Griis. & Gallagheil l2003l: Ide Griis et aLll2005h to de- 
termine the best match of the observed colours to the models and 
hence the ages of the open clusters: 

2 / ohsi — modeli \ 

^ =n — — ) 

where N is the number of colours (1-4) used in the analysis, obsi is 
the observed colour, models is the corresponding model colour, and 
(T,; is the uncertainty in the obsi colour. All ages with a fit of ^ N 
were considered good fits. For each colour set there was therefore 



a range of predicted ages for each OC. The age associated with the 
minimum is taken as the best-fit age. To determine the uncer- 
tainties in the predicted age we find the minimum and maximum 
ages within a A y^ defined to give a 68% confidence level (e.g., 
IPress et al]|l992h . Additionally, we add the model step size to the 
age uncertainty. We have also predicted the amount of extinction in 
each case. The amount of extinction applied to the model associated 
with the best-fit age is the best-fit E(B— V). The uncertainty in the 
predicted E(B— V) is determined from the minimum and maximum 
E(B— V) within the same Ax^ mentioned above. Additionally, we 
add the model step size to the E(B— V) uncertainty. 

In Table 1, we list the subset of our sample that has the U, B, 
and V magnitudes available. Column one is the name of the OC, 
column two is the published age, column three is the predicted 
age determined from the (U— B) and (B— V), column four is Mv, 
columns five and six are the reddening-corrected measured and zero 
reddened model (U— B) colours respectively, column seven is the 
difference of the measured and model (U— B), and columns eight, 
nine and ten are the same but for the (B— V) colour. It can be seen 
from Table 1 that for 74% of the OCs, the (U-B) differences are 
< 0.14, the assumed maximum measurement uncertainty in mea- 
sured colour, while 60% of the (B— V) differences are < 0.14. 



5 RESULTS 

5.1 Predicting Ages 

In Figure 2, we plot the predicted ages against the published ages, 
if only (U— B) and (B— V) are used to estimate the ages. We also 
provide a histogram showing the distribution of the differences be- 
tween the predicted and published ages, a plot of the predicted 
E(B— V)s against the published E(B— V)s, and a histogram show- 
ing the distribution of the differences between the predicted and 
published E(B— V)s. In this Figure, we used the solar abundance 
models. Plots of the predictions determined with the other colour 
combinations and at other assumed metallicities are similar so are 
not shown. We see from this figure that with the combination of 
(U— B) and (B— V) the ages of the OCs can be predicted with rea- 
sonable accuracy. 

We present Table 2 to discuss both the accuracy and preci- 
sion of various colours as age indicators. Column one is the colour 
combination, column two is the number of OCs in the sample with 
those colours available, and column three is the number of good 
fits for the solar, 0.2 solar and 2.5 solar models, where the number 
of good fits is the number of models in our grid that give ^ N. 
Column four is the percentage of recovered ages, which describes 
how often the predicted age uncertainties overlapped with the un- 
certainties in the measured ages when solar metallicity is assumed. 
The results from comparisons to models of 0.2 solar and 2.5 solar 
abundances are similar so are not shown. Column five lists the av- 
erage uncertainty in predicted log(t) in the negative direction for all 
the OCs, column six lists the average uncertainty in predicted log(t) 
in the positive direction for all the OCs, and column seven lists the 
average of the positive and negative uncertainties for all the OCs. 
It should be noted that the model ages range from a log(t) of 6 to 
10.3, therefore 2.15 is the maximum mean positive and negative 
uncertainty in log(t). Columns eight, nine, and ten are the percent- 
ages of predicted ages that are within a factor of 3 of the published 
ages for the solar abundance model, the 0.2 solar abundance model, 
and the 2.5 solar abundance model, respectively. This percentage is 
determined for the best-fit ages without consideration of the pre- 
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Table 1. Ages and Colour Differences 



Name 


Age^ 


Age^ 


Mv- 


(U-B) 


(U-B) 


diff. 


(B-V) 


(B-V) 


diff. 




(pub) 


(pred) 




measured 


model 




measured 


model 




Bochum 10 


8±2 


1 +31 
-0 


-6.52 


-1.05 


-0.67 


-0.38 


-0.24 


0.27 


-0.51 


NGC 6871 


9±2 


^-'-22 


-7.28 


-0.89 


-0.80 


-0.09 


-0.24 


0.17 


-0.41 


Stock 14 


10±2 


T , +96 
-'^-30 


-8.25 


-0.39 


-0.76 


0.37 


0.34 


0.17 


0.17 


King 12 


11±1 


, -,-^67 
^-'-12 


-4.90 


-0.88 


-0.72 


-0.16 


-0.20 


0.14 


-0.34 


IC 2581 


13±3 




-8.11 


-0.33 


-0.67 


0.34 


0.07 


0.05 


0.02 


NGC 3105 


21±3 


^+49 
^-3 


-6.95 


-0.69 


-0.67 


-0.02 


0.11 


0.06 


0.05 


NGC 6250 


22±5 


^^+157 
66 


-4.69 


-0.59 


-0.66 


0.07 


-0.16 


0.06 


-0.22 


NGC 6396 


30±8 


-20+2862 
3S_37 


-5.03 


-0.64 


-0.58 


-0.06 


-0.05 


0.07 


-0.12 


Trumpler 1 


30±6 


, , +2Q 
1 1 ^ 


-5.65 


-0.88 


-0.58 


-0.30 


-0.04 


0.07 


-0.11 


Lynga 6 


35±9 


30_2g 


-4.38 


-0.65 


-0.54 


-0.11 


0.01 


0.07 


-0.06 


Melotte 20 


43±18 


01 + 114 


-5.71 


-0.42 


-0.48 


0.06 


0.02 


0.07 


-0.05 


King 10 


45±11 




-5.34 


-0.58 


-0.47 


-0.11 


-0.19 


0.07 


-0.26 


NGC 129 


62±15 


n+20 


-6.03 


-0.64 


-0.40 


-0.24 


0.39 


0.09 


0.30 


NGC 6834 


65±18 


, T , -1-1 2fiQ 


-6.20 


-0.45 


-0.39 


-0.06 


-0.10 


0.10 


-0.20 


NGC 6405 


71±21 


, ^ , -1-1 4-fi 


-3.65 


-0.27 


-0.38 


0.11 


-0.03 


0.11 


-0.14 


Collinder 394 


74±6 


0-7 -1-178 


-4.74 


-0.26 


-0.38 


0.12 


0.23 


0.11 


0.12 


NGC 6025 


74±22 


^^4-191 

82^121 


-4.38 


-0.46 


-0.38 


-0.08 


-0.12 


0.11 


-0.23 


NGC 5662 


77±20 


222 


-4.38 


-0.24 


-0.37 


0.13 


0.30 


0.11 


0.19 


NGC 6087 


78±19 


, .__l_309 
^^^-137 


-5.09 


-0.21 


-0.37 


0.16 


0.18 


0.12 


0.06 


NGC 1778 


129±29 




-4.16 


-0.27 


-0.25 


-0.02 


0.09 


0.14 


-0.05 


NGC 1647 


130±25 


1 1 1 -1-197 
^*^-105 


-2.35 


-0.35 


-0.25 


-0.10 


0.07 


0.14 


-0.07 


NGC 5316 


166±33 


1 rt-l-12681 
^*^-113 


-4.91 


0.13 


-0.19 


0.32 


0.67 


0.15 


0.52 


NGC 744 


184±49 


^-'"-170 


-5.17 


-0.14 


-0.16 


0.02 


0.06 


0.15 


-0.09 


Melotte 105 


224±53 




-4.27 


-0.11 


-0.10 


-0.01 


0.07 


0.16 


-0.09 


NGC 3532 


262±46 


A J —1— 3 3 

24315^^ 


-5.50 


-0.03 


-0.04 


0.01 


0.22 


0.18 


0.04 


NGC 2548 


364±102 


^r.^-l-411 


-3.23 


0.12 


0.07 


0.05 


0.26 


0.23 


0.03 


Melotte 111 


COO 1 


384_2i7 


-1.95 


0.12 


0.16 


-0.04 


0.20 


0.33 


-0.13 


NGC 2527 


619±163 


T 1 '7 -|— 783 
•"'-204 


-1.90 


0.09 


0.18 


-0.09 


0.30 


0.39 


-0.09 


NGC 2266 


lib±ll 


/IQO+1302 
498 _ 348 


-3.83 


0.40 


0.18 


0.22 


0.62 


0.46 


0.16 


Berkeley 2 


794dil 


999 _ 953 


-4.18 


0.14 


0.17 


-0.03 


0.46 


0.49 


-0.03 






ToS+3562 


-J.O / 




n 1 "7 
U. 1 / 


U.Uj 


Vi.jy 


n ^ 1 
U.J 1 


o oc 


Berkeley 69 


867±48 


506118594 


-3.16 


0.21 


0.17 


0.04 


0.38 


0.52 


-0.14 


NGC 2477 


875±238 


9cr,+ 10350 
^-»J_180 


-5.70 


0.17 


0.17 


0.00 


0.49 


0.53 


-0.04 


NGC 2192 


1072±48 




-3.53 


0.10 


0.17 


-0.07 


0.49 


0.62 


-0.13 


NGC 2660 


1351±291 




-4.08 


0.34 


0.22 


0.12 


0.61 


0.72 


-0.11 


NGC 1798 


1421±16 


47+9653 
^'-44 


-4.86 


0.15 


0.25 


-0.10 


0.82 


0.77 


0.05 


NGC 2506 


1648±485 


2^9+2941 
^J^_175 


-4.31 


0.22 


0.28 


-0.06 


0.54 


0.82 


-0.28 


NGC 7044 


1824±361 


lQQQQ+102 

'""-19810 


-3.93 


0.48 


0.26 


0.22 


0.58 


0.81 


-0.23 


Berkeley 32 


3477±698 




-2.97 


0.28 


0.36 


-0.08 


0.78 


0.91 


-0.13 


NGC 6253 


3949±1086 


6800ti?399 


-2.68 


0.34 


0.38 


-0.04 


0.81 


0.92 


-0.11 


NGC 2682 


4093±958 


,3QQ+8200 
lJUU_j^294 


-3.16 


0.28 


0.38 


-0.10 


0.78 


0.92 


-0.14 


NGC 6791 


7850±2026 


1QQQ9+102 
'""-19805 


-4.14 


0.82 


0.46 


0.36 


1.02 


0.97 


0.05 



1 Ages in Myr from lPaunzen & Netopij 12006) 
^ Ages in Myr predicted from (U— B) and (B— V) 



dieted or measured age uncertainties. Unfortunately, the only sam- 
ples of OCs large enough to draw any significant interpretations are 
the (U-B), (B-V) and (U-B) and (B-V). Results from the other 
combinations of colours are merely suggestive. 

We see from Table 2 that for several colour combinations the 
uncertainties in predicted age are relatively large. This is because 
we have assumed the maximum uncertainty in measured colours 
for the OCs. Using the actual uncertainties in integrated colours for 
this sample would not change the predicted ages but would result 
in better age constraints. Unfortunately the uncertainties in inte- 



grated colours of each of these OCs are not available. However, 
the relative predicted age uncertainties provide a good measure of 
the precision afforded by each colour and colour combination for 
age-dating. Looking at the mean uncertainties (column 7 in Table 
2) we see that the precision is improved when the (U-B) colour is 
included. 

We see from this table that by using the (U-B) and (B-V) 
colours up to 83% of the predictions are within a factor of 3 of the 
puWished values. For the best-fit (U-B), (B-V), and (U-B) and 
(B-V) samples, we have determined the differences between the 
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Figure 2. Top Left: Model predicted age versus the published age of the OCs in the WEBDA sample, as determined by the (U— B) vs. (B— V) colours. The 
solid black line is y=x, the inner dashed lines represent a factor of 3 difference in age, and the outer long-dashed lines represent a factor of 10 difference in 
age. Top Right: Histogram showing the distribution of the difference between the predicted and published ages. The inner and outer dashed lines represent a 
factor of 3 and 10 change in age respectively. Bottom Left: Model predicted extinction vs. published extinction. The solid black line is y=x, the inner dashed 
lines represent a change of 0.2 mag, and the outer long-dashed line represents a change of 0.5 mag. Bottom Right: Histogram showing the distribution of the 
difference between the predicted and published E(B— V). The inner and outer dashed lines represent a change of 0.2 and 0.5 mag respectively. 



best-fit predicted and measured log ages. We find that 68% of the 
differences are less than 0.77 dex, 0.86 dex, and 0.44 dex for the 
(U-B), (B-V), and (U-B) and (B-V) samples respectively. This 
indicates that (U-B) alone gives more accurate ages than (B-V) 
alone. We also see that combining (U-B) with (B-V) further im- 
proves the accuracy of the predicted ages. 

Table 2 shows the importance of the (U— B) colour in measur- 
ing the age of a cluster. To the extent that our sample size allows 
such comparisons, combinations of colours that include (U-B) 
tend to have more accurate and more precise age predictions. The 
fact that single colours do not do a good job of constraining the ages 
is likely because there are three degenerate unknowns, age, redden- 
ing, and metallicity, and two filters are not enough t o break these 
degen eracies. This is consistent with results from I Anders et al.l 
( I2OO4I) . However, the (U-B) colour alone is a better predictor of 
age than any combination of colours that do not include (U-B), 
regardless of the assumptions about metallicity, consistent with re- 



sults from iGil de Paz & Madorg | |2002|) and lAnders etaU ( l2004t) . 
The (U-B) colour reflects the contributions from the young and 
intermediate age stars. Tracing the younger stars in a population is 
critical to constrain the age. Colours that do not include U only fol- 
low the slower-evolving stars and are not as sensitive to changes in 
age. 



5.2 Extinction and Metallicity 

The extinction is predicted well by the same colours that predict 
the age well. For example, using the (U— B) and (B-V) colours 
we find that 62% of the E(B— V) estimates are within 0.2 mag of 
the published values. The colour that is the most sensitive to age, 
(U— B), is once again critical. 

We get similar numbers of good fits to the cluster ages whether 
we assume 0.2 solar, solar, or 2.5 solar abundances (see Table 2). 
This implies that, at least within the metallicity range of our sample 
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Table 2. Prediction Quality 
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Figure 3. Published My plotted against the pubhshed ages of the open 
clusters in our sample. The curve app roximates the lowest luminosity limit 
curve in lCerviiio & Luridianal ^20041) . The filled black squares ai'e the OCs 
fainter than the curve while the filled black triangles are brighter. 

OCs, knowing the metallicity is not critical for getting a reasonable 
fit to the cluster age. Note that we are not fitting for the metallicities 
of these clusters, but instead are testing to see how well we can 
recover the correct ages and extinctions with standard solar and 
near solar metallicity models. 

5.3 The Effect of Hq 

The SB99 model colours do not include nebular Ha emission. This 
does not impact our analysis of the integrated colours of the OCs 
because the stellar components are resolved and no Ha is included 
in either the models or the data. However, this should be considered 
when models are compared to young unresolved clusters in other 
galaxies. 

For clusters younger than ~10 Myr, the Ha emission con- 



tributes significantly to the Johnson R band and minimally to the 
Johnson V band. From SB99, it is seen that a 1 Myr model has an 
EW(Ha) of ~2500 A, a 5 Myr model has an EW(Ha) of ~400 A, 
and by 10 Myr the EW(Ha) is ~10 A. For comparisons to young 
unresolved clusters. Ha emission should be included in the model 
SEDs. For a discussion about adding Ha emission to m odels to 
compare to the clusters in Arp 285, see ISmith et al.l ( l2008h . 



5.4 The Effect of Stochastic SampHng of the IMF 

It has been shown that traditional synthesis models usually re- 
turn the mean value of the integrated stellar population distribu - 
tion (e.g. ICervino & Luridian3l2004 ITuridiana & Cervindl2007h . 
While this is correct on average, it is not necessarily correct in in- 
dividual cases because the mean may not be representative of the 
actual values. Galactic open clusters are typically low mass so will 
be subject to the effects of stochastic sampling of the IMF (e.g. 
ICervino & Valls-Gabau j|20 03). 

Cervino & Luridianal j2004) define the lowest luminosity limit 
as a criterion that must be met in order to compare an observed clus- 
ter to a model cluster. The total luminosity of the observed cluster 
must be larger than the indi vidual contribution of any of the stars 
in the model. According to ICervino & Luridianal ( |2004|) . clusters 
brighter than this limit may or may not have a well-sampled IMF, 
but clusters fainter than this limit will be misrepresented by synthe- 
sis models. 

Figure 3 is a plot of the published M\/ plotted against the pub- 
lished ages of the OCs in our sa mple. The curve approximat es the 
lowest luminosity limit curve in ICerviiio & Luridianal ( |2004|) . The 
filled squares are the OCs fainter than the curve, while the filled 
triangles are OCs that are brighter. There are 15 OCs fainter and 
27 OCs brighter than the curve. From this figure it can be seen that 
several of the OCs in our sample are brighter than the lowest lumi- 
nosity limit, but the majority are fainter. This suggests that many 
clusters in our sample (those fainter) should show these stochastic 
effects, while several might not (those brighter). 

Figure 4 shows the distributions of the differences between 
the published and the model (U-B) colours and the published and 
model (B— V) colours for the OCs in our sample. The measured 
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Figure 4. Top: Distribution of tlie difference between the publislied 
(U— B)p„i, and the model (U— B)^^^ for the open clusters in our sample. 
The model colours assume the published ages. The colours are reddening- 
corrected. The histograms are divided into two groups; fainter than and 
brighter than the lowest luminosity hmit curve, as shown in Figure 3. Bot- 
tom: Same as above but for the (B— V) colours. 



Figure 5. Top: Difference between the published and model (U— B) colours 
against the published ages of the clusters. Bottom: Difference between the 
published and model (B— V) colours against the published ages of the clus- 
ters. The two panels describe the same two groups as above. The dashed 
lines show zero colour difference. 



colours are reddening-corrected. The model colours assume the 
published ages. To make comparisons, only OCs with both inte- 
grated (U— B) and (B— V) colours are included in these distribu- 
tions. The histograms are divided into 2 groups, those brighter and 
those fainter than the lowest luminosity curve, as in Figure 3. 

For the (U— B) colour, a Kolmogorov-Smimov (K-S) test can- 
not rule out that the brighter and fainter samples originated from 
the same parent distributions. The rms for the differences in (U— B) 
are 0.16 and 0.17 for the brighter and fainter samples respectively. 
However, for (B— V), a K-S test suggests a probability of only 
0.012 that the two samples (brighter and fainter) originate from 
the same parent distribution. The rms values are 0.15 and 0.25 for 
the (B— V) distributions respectively. The scatter appears slightly 
larger for the fainter samples, with the (B— V) difference being 
slightly skewed to negative numbers (see Figure 4). 

Figure 5 shows the differences between the published and 
model (U— B) and (B— V) colours plotted against the published 
ages using the same sample as in Figure 4. The measured colours 



are reddening-corrected. For each colour difference the OCs are 
divided into the same bins as described above. From Figure 5 it 
can be seen that the scatter in the (U— B) colour differences do 
not appear to be correlated with age. The fainter (B— V) differ- 
ence sample appears to have a slight dependence on age in that 
the published (B— V) colours are bluer th an the model co lours for 
the younger OCs, consistent with results in lBruzuail ( |2002|) (see be- 
low). Given this, the K-S results, and the larger rms, we might be 
seeing stochastic sampling effects in the differences between the 
measured and model (B— V) colours. 

In Figure 6, we plot the published reddening-corrected (U— B) 
colour against the published reddening-corrected (B— V) colour for 
our OCs. This figure is divided into the same two groups described 
above. Also on this figure are our SB99 models with E(B— V) of 0.0 
mag (solid curve) and 0.5 mag (dashed curve). This figure shows 
that the scatter of these colours around the models is not a function 
of Mv'. Note that relatively young clusters (that is, blue clusters) 
tend to lie above and to the left of the zero extinction model curve. 
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Figure 6. Published (U-B) against published (B-V) for the OCs. The fig- 
ures are divided into 2 groups as described above. The curves are our SB99 
models with E(B— V) of 0.0 mag (solid curve) and 0.5 mag (dashed curve). 



Figures 4 and 5 show that this is predominately due to the fainter 
clusters being bluer in (B— V) than the models. 

Our results can be compared to population synthesi s studies 
that ex plici tly investigate stochastic e ffects. For example, ISruzuall 
( I2OO2I) and lBruzual & CharlotI j2003l) have run Monte Carlo simu- 
lations of low mass star clusters, to investigate how stochastic ef- 
fects can cause observed colours to differ from models with full- 
sampled IMFs. In their simulations, young model clusters with 
masses of 10^ Mq also lie a bove models w ith fully-sampled IMFs 
on a (U-B) vs. (B-V) plot |l Bruzualll2002i) . consistent with what 
we see in Figure 6. Thus at least some of the observed scatter may 
be due to stochastic effects. 

From these figures and tables we see suggestive evidence of 
stochastic sampling of the cluster IMF in this sample of OCs. It 
is possible that the uncertainties in the measured colours and/or 
ages are large enough (0.14 mag and 19% respectively) that the ef- 
fects of stochastic sampl i ng are washed out. However, figure 1 in 
ICervino & Valls-Gabau j ( l2003l) shows that the expected scatter in 
the integrated (B-V) of a 10 Myr cluster with 10"^ or less stars 
is much larger than our assumed colour uncertainties. Most of the 
OCs in our sample are much older than 10 Myr. It is possible that 
stochastic sampling effects are smaller in (B— V) for older clus- 
ters. After ~ 10 Myr the most luminous stars will have evolved 
off of the main sequence and will co ntribute more strongl y to the 
red end of the integrated cluster SED. Ide Griis et aL U2OO5I) " remark 
that stochastic sampling effects affect broadband photometry to a 
smaller extent than spectroscopy. In a follow up paper (Hancock 
et al., in prep.) we will explore a more detailed analysis of the 
stellar content of these resolved OCs to further investigate the ex- 
tent to which stochastic sampling of the IMF affects the integrated 
broadband colours. If stochastic sampling effects are observation- 
ally characterized in this sample of OCs, they will provide a good 
bench mark for testing future population synthesis models. 

5.5 The Effect of Cluster Dissolution 

Current stellar population synthesis models do not take into ac- 
count the effects of cluster dissolution. Mass segregation causes 



the massive stars to concentrate towards the centre of a clus- 
ter while the low-mass stars tend to populate the outer regions. 
A consequence of mass segregation is that low-mass stars tend 
to be preferentially ejected from star clusters. The mass func- 
tion (MF) therefore changes drastically during the evolution of 
a cl uster, as a consequence, bec oming less steep than the IMF 
(e.g. iBaumgardt & Makind |2003|) . The preferential loss of low- 
mass stars will change the integrated colours of a cluster (e.g. 
iLamers. Anders. & de G riis 2006). Such a time evolution of the in- 
tegrated colours resulting from cluster dissolution can be difficult 
to disentangle from stellar evolutionary effects. 

The colour change of a star cluster due to dissolution de- 
pends upon its age, its mass, and its dissolution time scale 
(tdia), with the effect being larger for later ages and longer l^ia ■ 
iLamers. Anders. & de GriisI ( l2006h show that, for a 10^ M© clus- 
ter, the (B-V), (V— I) and (V— K) colours are not significantly af- 
fected by the preferential loss of low-mass stars during the first 80% 
of their tdis, with A(V— 1)= 0.03 and 0.1 mag for a Xdis of 1 Gyr 
and 10 Gyr respectively. However, in the last 20% of a cluster's life, 
the colour effect can be quite large, up to A(V— 1)= 0.3 mag. 

Our sample clusters likely have masses between 10'' and 
10^ M(7). This implie s idis between ~100 Myr and 2 Gyr 
jLamers & Gielesll2006l) . Given that 45% of our sample have ages 
less than 80 Myr and 86% have ages less than 1.6 Gyr (see Table 1), 
it is likely that most of our clusters are not in the last 20% of their 
lifetimes, but we can not rule out the possibility that some might 
be near the end of their lifetimes. Thus the effect of dissolution 
on their colours is likely small, typically smaller than our assumed 
uncertainties of 0.14 magnitudes. 

Using the stellar data available on WEBDA for a sample of 
our OCs covering a representative range of distances, we found that 
~ 85% contained stars <0.5 Mq. To determine the effects that 
low mass stars have on the model integrated colours we compared 
Kroupa models with lower mass cutoff of 0. 1 Mq with those with 
0.5 Mq . We found no significant differences in the (U— B), (B— V), 
(V-R), or (V-I) colours. 



6 SUMMARY 

We present an empirical assessment of the use of broadband colours 
as age indicators for unresolved extragalactic clusters and investi- 
gate stochastic sampling effects on integr ated colours. T he popu- 
lation synthesis code Starburst99 tLeitherer et and four 
optical colours were used to estimate how well we can recover the 
ages of 62 well-studied Galactic open clusters with published ages. 
We conclude the following: 

1) Galactic open clusters can be used for testing the inte- 
grated properties from population synthesis and serve as reasonable 
benchmarks for future assessments of age-dating methods. 

2) The (U— B) colour is critical. Only colour combinations that 
included (U— B) resulted in good age and extinction predictions, 
consistent with previous results. 

3) Only with (U-B) included were the predicted age uncer- 
tainties reasonably constrained. 

4) Changes of a factor of ~2 in assumed metal abundance 
do not result in significantly different predictions of cluster age. 
This indicates that the uncertainties in predicting the age of a clus- 
ter resulting from the age-reddening degeneracy dominate over the 
other sources of degeneracy in these optical bands, over our age 
range and metallicity range. Another possibility is that uncertain- 
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ties in the measured colours (and hence the predicted ages) dwarf 
the metallicity-age degeneracy. 

5) A minimization and a Ax^ defined to give 68% confi- 
dence levels provide reliable age estimates, and more importantly, 
reliable age uncertainties. The difference in the photometric ages 
and the ages derived from the HR-diagrams of our selected cluster 
sample with both the (U— B) and (B— V) colours are smaller than 
0.5 dex for 79% of the clusters, with a maximum difference of 1.5 
dex. 

6) It is likely that the uncertainties in the measured colours 
and/or ages are large enough (0.14 mag and ~ 20% respectively) 
that the effects of stochastic sampling are washed out. A more de- 
tailed analysis of the stellar content of these resolved OCs will be 
required to further investigate the extent to which stochastic sam- 
pling of the IMF affects the integrated broadband colours. 

7) If stochastic sampling effects are observationally character- 
ized in this sample of OCs, they will provide a good bench mark 
for testing future population synthesis models. 
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